Finite-size effects in the self-organized critical forest-fire model 
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We study finite-size effects in the seif-organized criticai forest-fire modei by numericaliy evaluating 
the tree density and the fire size distribution. The results show that this model does not display the 
finite-size scaling seen in conventional critical systems. Rather, the system is composed of relatively 
homogeneous patches of different tree densities, leading to two qualitatively different types of fires: 
those that span an entire patch and those that don't. As the system size becomes smaller, the 
system contains less patches, and finally becomes homogeneous, with large density fluctuations in 
time. 



I. INTRODUCTION 

During the past years, systems which exhibit self- 
organized criticality (SOC) have attracted much atten- 
tion, since they might explain part of the abundance 
of fractal structures in nature |^]. Their common fea- 
tures are slow driving or energy input and rare dissipa- 
tion events which are instantaneous on the time scale of 
driving. In the stationary state, the size distribution of 
dissipation events obeys a power law, irrespective of ini- 
tial conditions and without the need to fine-tune param- 
eters. Examples for such systems are the sandpile model 
0, the self-organized critical forest fire model ||,^, the 
earthquake model by Olami, Feder, and Christensen [Q, 
and the Bak-Sneppen evolution model Numerical 
as well as analytical studies of those systems are usually 
based on the assumption that their critical behaviour can 
be described in similar terms as that of equilibrium crit- 
ical systems. This assumption is given a basis in 
where it is suggested that SOC systems can be mapped 
on conventional critical systems by interchanging control 
and order parameters. Thus, the Bak-Sneppen model 
can be mapped on a depinning problem However, it 
has been shown in |^] that the mapping suggested in ||] 
for the SOC forest-fire model does not generate a system 
with a conventional critical point. Instead, the phase 
transition shows hysteresis effects and is discontinuous 
when approached from above. Other unconventional fea- 
tures have also been seen in the SOC forest-fire model, 
like the existence of more than one diverging length scale 
1^,^, the absence of a spanning cluster immediately be- 
yond the critical point and the dependence of the 
large-scale behaviour on details of the model rules |P,[lO| . 
(For a review on the SOC forest-fire model, see [[13|.) 
Other SOC systems show also unconventional scaling be- 
haviour. Thus, in the two-dimensional abelian sandpile 
model finite-size scaling is violated and the critical 
exponents for the earthquake model by Olami, Feder, and 
Christensen Q appear to depend continuously on the pa- 
rameters. There is substantial need to better understand 



the nature of the scaling behaviour of those systems. 

It is the purpose of this paper to shed some light on 
the unconventional critical behaviour of the SOC forest- 
fire model by studying its finite-size effects. We choose a 
version of the model which is identical to the SOC forest- 
fire model for system sizes much larger than the correla- 
tion length, and we discuss the changes that occur in the 
model as the system size is decreased below the correla- 
tion length. We find that instead of displaying finite-size 
scaling, small systems undergo a rearrangement from a 
structure with patches of different density to a more ho- 
mogeneous structure with large density fluctuations in 
time. We find also that, contrary to conventional critical 
systems, small systems and small parts of large systems 
differ in the probability distribution for the density and 
in the fire-size distribution. We suggest that these re- 
sults can be explained by the fact that the system has 
two qualitatively different types of fires. 

The outline of this paper is as follows: In section II, 
we define the model that we used for studying finite-size 
effects, and discuss briefly known results. Section III 
shows computer simulation results for the fire-size distri- 
bution and the tree density as the system size changes 
from values larger than the correlation length to values 
much smaller than it. In the conclusion, we summarize 
and discuss our findings. 



II. THE MODEL 

The version of the SOC forest-fire model studied in this 
paper is defined on a square lattice with sites. Each 
site is either occupied ("tree") or empty ("no tree"). At 
each time step, the system is updated according to the 
following rules: (i) "Burning" : A site in the system is 
chosen at random ("struck by lightning"). If the site is 
occupied, the whole cluster of occupied sites connected 
to this site (by nearest-neighbour coupling) is removed 
from the system ("burnt"), i.e., the occupied sites of that 
cluster turn to empty sites. If the chosen site is empty, 
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nothing happens, (ii) "Tree growth" : We select randomly 
So = pL^ sites from the system and occupy those that are 
empty (possibly also including sites which have become 
empty due to the removal of the cluster). These sites are 
selected one after another, allowing for the same site be- 
ing selected more than once during the same filling step. 
In principle, sq can therefore be larger than L^, however, 
in our simulations we chose usually values smaller than 

For fixed so and very large system size L, these rules 
are equivalent to having a lightning probability f — 1 / L"^ 
per site and time step, and a tree growth probability p, 
and the model is identical to the original SOC forest-fire 
model H]. Because of this equivalence, which was first 
pointed out by Grassberger most numerical studies 
of the SOC forest-fire model up to now were performed 
using the above rules, which allow for fast and efficient 
computer simulations. With the above rules, finite-size 
effects can also be studied very efficiently, as was sug- 
gested in However, one has to keep in mind that the 
results are somewhat different from those for the original 
model. While in the original model lightning can strike 
the system between the growth of any two trees, it can 
strike the system in the present model only after growth 
step (ii) is finished. This leads to density peaks in Fig- 
ures ^ and |l^ below that are not present in the original 
model. However, our main conclusions are not affected 
by the particular choice of the dynamical rules, as will 
be discussed further below. 

Let us first summarize shortly the major numerical re- 
sults for the case so ^ L^, as reported in the literature on 
the SOC forest-fire model In this hmit, only 

a small number of trees grow at each time step (com- 
pared to the total number of trees). After a transient 
time, a stationary state is reached where the tree density 
has only small fluctuations around some average value 
p(so) that does not depend on L. Throughout this pa- 
per, we study only stationary states and do not evaluate 
the initial transient behaviour. Since the mean number of 
trees s burnt during a fire must be identical to the mean 
number of trees growing between two fires, we have the 
relation 



s = so(l - p)/p. 



(1) 



The leading finite-size corrections to this equation are of 
order sq/L'^ and can be neglected in the case Sq <^ L? 
which we are considering in this paragraph. As So in- 
creases, the mean fire size increases also, and we approach 
the critical point of the SOC forest-fire model, where the 
mean tree density is given by pc — 0.41. The correlation 
length ^ is a measure for the radius of the largest tree 
cluster and is related to sq via ^ ~ Sq, with « 0.58 in 
d = 2 dimensions. The size distribution of tree clusters 
near the critical point is well described by the scaling 
form 



n(s) 



T(s/s„ 



(2) 



with a cutoff function C that is constant for small argu- 
ments and decays exponentially fast when the argument 
is considerably larger than 1. The cutoff cluster size Smax 
is related to the correlation length ^ via Smax S,'^ , with 
fi being the fractal dimension of tree clusters, which is 
found to be 1.95 or 1.96 |Jll|. The value of the 
exponent t is approximately 2.14. The relation between 
Smax and So is s^ax ^ Sq' ■^it^^ A = ~ 1.15 

All these numerical findings agree well with conven- 
tional scaling assumptions based on a single diverging 
length scale. Analytical studies of the model, such as 
mean-field theories jr^-p^ and rcnormalization group 
calculations ||2^,|2l[ are also based on conventional scaling 
assumptions. Therefore, the violation of finite-size scal- 
ing described in the following might appear surprising 
to many readers. However, one must keep in mind that 
the simulation data do not cover much more than one 
decade in the correlation length ^. The observed scaling 
behaviour eq. (||), together with the measured values of 
the critical exponents, do not necessarily indicate an ex- 
act asymptotic scaling form, but may simply be a good 
approximation to more complicated scaling, which works 
well for the system sizes and parameter values stud- 
ied in simulations. A similar phenomenon is known for 
the sandpile model, where good scaling collapses for the 
avalanche size distribution could be achieved in 23 
and older papers, although it has been recently shown 
iQ that finite-size scaling is violated and that the sim- 
ple scaling ansatz used for the data collapse is incorrect. 

Figure ^ shows a snapshot of a system with a tree 
density p just below pc- One can distinguish regions of 
different densities with a rather homogeneous tree dis- 
tribution within a region. These regions are obviously 
created by a fire that burns down a cluster of high tree 
density. After the fire, a burnt region is almost empty 
and becomes slowly filled with trees according to the law 
p — p{l — p). We call these regions of homogeneous tree 
density "patches" , as we did in . If lightning strikes 
a tree in a patch of low density, it usually burns down a 
small tree cluster. If it strikes a patch of a density larger 
than the percolation threshold, it burns down a tree clus- 
ter as large as the patch itself. This observation indicates 
that there are two qualitatively different types of fires in 
the system: those that span an entire patch, and those 
that destroy a small percolation cluster within a patch 
of a tree density below the percolation threshold. As we 
will see below, this gives rise to the unusual finite-size 
properties of the model. 

If the correlation length ^ is of the same order as or 
larger than L, the behaviour sketched above is modi- 
fied due to finite-size effects. For not too small values of 
sq/L^, the tree density increases by a noticeable amount 
between two fires, leading to large density fluctuations 
and to fires that span the entire system. If the SOC 
forest-fire model showed conventional critical behaviour, 
there would be a single diverging length scale, namely 
the correlation length ^, which would be related to f /p 
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or, equivalently, to sq, via ^ ~ {f/p)~'^ or ^ ~ Sq. Finite- 
size effects would then manifest themselves in a scaling 
form 

n{s)^s-^C{s/L''), (3) 

for the size distribution of tree clusters. Furthermore, 
on scales smaller than L and ^, all measured quanti- 
ties should be indistinguishable from those measured in 
a small section of an infinitely large critical system. 

The following section presents simulation results that 
show that none of these finite-size scaling assumptions 
is satisfied for the SOC forest-fire model. In fact, the 
invalidity of the assumption of a single diverging length 
scale has already been shown in |^ . The invalidity of the 
second assumption that measurements in small systems 
and in small sections of large systems should give identi- 
cal results, can be understood by considering for instance 
the mean time interval between two fires. In a small sub- 
system of linear size / of a large system with a correlation 
length ^ ^ I, this is given by (p(l — pc))~^. Just before 
fire reaches the subsystem, its tree density is far above 
the percolation threshold, and the spanning cluster of the 
subsystem is part of a large tree cluster that extends far 
beyond the limits of the subsystem. Lightning usually 
strikes this large cluster outside the subsystem, the time 
interval between two lightning strokes within the subsys- 
tem being L'^/P, which diverges as L diverges. In con- 
trast, fire cannot enter a small system from outside, but 
the tree density of a small system increases until light- 
ning strikes a tree within the system. According to our 
rules, time is measured in units of the mean time inter- 
val between two lightning strokes. On this time scale, 
the time between two fires within a small system of lin- 
ear size I is finite. In contrast, the time interval between 
two fires within a subsystem of size I of a much larger 
system is vanishingly small compared to the time inter- 
val between two lightning strokes within the subsystem. 
All these arguments are backed up and complemented by 
the numerical results reported in the following. 

III. RESULTS OF COMPUTER SIMULATIONS 

In this section we will present and explain data ob- 
tained from about 300 runs of the model for various val- 
ues of So and L. Since many runs of the simulation were 
necessary, we chose a cluster of workstations rather than 
a " supercomputer" . The system size L varied between 
10 and 2000 in these runs. We found that as finite- 
size effects become more important, the system shows a 
transition between two qualitatively different types of be- 
haviour which we call critical behaviour and percolation- 
like behaviour. The critical behaviour is is characterized 
by a good scaling collapse of the fire size distribution 
and by large spatial variations in the local tree density. 
The percolation-like behaviour is characterized by large 
temporal fluctuations in the global tree density, with a 



rather homogeneous tree distribution within the system 
for any given time. Snapshots of the system therefore 
resemble percolation systems where each site is occupied 
by a tree with a probability p. (For an introduction to 
percolation theory, see e.g. pj|.) The following three 
subsections show how this transition manifests itself in 
the mean tree density, the fire size distribution, and the 
probability distribution for the tree density. 



A. Lines of constant tree density 

First, we measured the mean tree density 

T 

p^{l/T)J2pit) (4) 

t=i 

in the system, averaged over a large number of T itera- 
tions, for various values of L and sq. The density p{t) was 
always evaluated after the refilling step (ii) . Compared to 
a model where trees grow at a rate p = p(l — p), the den- 
sity values in our model are somewhat larger when sq / 
is not very small. If, for instance, the density is increased 
from p — Ap to p during the refilling step (ii) , the value p 
enters the above sum, while an evaluation based on a con- 
stant growth rate would give 1 — Ap/ ln(l + Ap/ (1 — p)) 
instead of p. 

It it not obvious what the relation between sq, and 
p should be if we want to deduce it from an analogy with 
equilibrium critical systems. We have already mentioned 
that the temporal fluctuations in p{t) become larger as 
the ratio Sq/L'^ increases. Similarly, temporal fluctua- 
tions increase in a critical equilibrium system when the 
system size becomes smaller. One might therefore expect 
that decreasing L at fixed so should drive the system to- 
ward the critical point, where p — pc — 0.41. However, 
we have argued in the previous section that a given site 
burns down more often in a large system without finite- 
size effects than in a smaller system with finite-size ef- 
fects that has the same value of sq. From this, it follows 
that the mean tree density increases with decreasing L, 
when So is fixed. In the limit so 3> it must go to 
one. From this point of view, a system with sufficiently 
large finite-size effects should rather be compared to an 
equilibrium system in the ordered phase, for instance to 
a percolation system beyond the percolation threshold. 
The analogue of sq in a- percolation system is then the 
mean size of the cluster that a given site belongs to, and 
it is proportional to beyond the percolation threshold 
p^ . Indeed, we find that the mean size of fires becomes 
proportional to when finite-size effects are strong (see 
below). However, there is nevertheless a fundamental 
difference between a percolation system beyond the per- 
colation threshold and our system with a density above 
Pc'- In our model, a tree cluster that spans the system and 
has a size proportional to occurs only rarely when the 
mean density is only slightly above the critical density. 



3 



while a percolation system above the percolation thresh- 
old has always a system spanning tree cluster. 

In Figure lines of constant mean tree density are 
plotted on a double logarithmic scale in sq and L^. This 
figure shows that there are two qualitatively different re- 
gions in the sq vs. plane, with a transition region 
between them. First, there is the region where there are 
no finite-size effects, sq ^ L^. In this region, the size 
of the tree clusters is much smaller than the system size, 
and the global density fluctuations are small. The mean 
tree density of such a system is smaller than pc « 0.41. 
For systems with small density fluctuations the proba- 
bility that a given empty site is filled with a tree during 
one time step, is given by sq/L'^, and the probability 
that a given tree is burnt by a fire is so(l — p)/(pL'^). 
Since both probabilities decrease as with increasing 
L, dynamics of larger systems are slower than those of 
smaller systems. Apart from this change of the charac- 
teristic time scales, the local dynamics is independent of 
L, and consequently correlation functions and cluster size 
distributions are the same for systems of different sizes 
(provided that ^ sq). This leads to the horizontal 
slope in the large L regime of the curves for p < pc m 
Figure ||. 

The curves show deviations from the horizontal be- 
haviour when L becomes as small as or smaller than the 
correlation length f ~ Sq. These deviation occur in the 
transition region where finite-size effects begin to become 
noticeable. The dynamics changes from fires that are not 
affected by the finite system size to a fire size distribu- 
tion that includes sometimes events of the order of the 
system size. Fires of such a large size destroy the patchy 
structure of the forest described earlier, and cause a more 
random tree distribution. 

In the second region, lines of constant mean density are 
curves of constant sq/L^. This feature can best be under- 
stood when considering the parameter range where sq is 
of the order of or above. In this range, the mean tree 
density is much larger than pc, and the system contains 
a spanning cluster after each filling. During the "burn- 
ing" step, this cluster is removed with a finite probability, 
leading to a large change in density in the system. When 
a "finite" (i.e. not spanning) cluster is removed during 
the "burning" step, the overall density hardly changes 
for large L. The time series of the density is therefore 
determined almost completely by the filling events and 
by the large burning events. Since the filling events fill 
a large fraction of empty sites, and since large burning 
events burn a large fraction of trees, the tree distribution 
within the system is rather homogeneous. A snapshot 
of the system at a given time looks therefore similar to 
a percolation system with a density p{t). From percola- 
tion theory we know that for a given density the fraction 
of trees sitting in the spanning cluster is independent 
of L, and consequently curves of constant large density 
are curves of constant sq / in Figure |[ Even curves for 
smaller s^jl?, which correspond to densities only slightly 
above pc show for sufficiently large L an asymptotic be- 



haviour So/Li^ = const, with the constant vanishing at 
Pc- The reason is again that finite fires do not reduce the 
density of an infinitely large system, and that the system 
spanning fires reduce the density by an amount that does 
not depend on L, but only on the density itself. 

Finally, the critical curve (for the density p — pc — 
0.41) is obtained from the condition 

i w ^ - Sq, with V ~ 0.58 . 

This curve is the separatrix between the two regions de- 
scribed above and is indicated in Figure ^ by the bold 
line. As L decreases, more and more curves merge with 
the separatrix when their correlation length becomes 
comparable to the system size. 

B. The fire size distribution 

Since lightning strikes each tree with the same prob- 
ability, the size distribution of fires is proportional to 
sn(s), with n{s) being the size distribution of tree clus- 
ters. As mentioned above, conventional scaling would 
imply a form sn[s) ~ s^~'^C{s / s'^'^) for the fire size dis- 
tribution if the correlation length ^ ~ Sq is smaller than 
the system size, and a finite-size scaling form sn{s) ~ 
s-^~'^C{s/ Lf^) in the opposite case. In both cases, one 
would obtain a scaling collapse of the curves for different 
So or L. 

Figure ^ shows the fire size distribution for parameters 
such that ^ < L. While not perfect, the data collapse is 
good and would not impose the conclusion that simple 
scaling is violated. The bump near the end of the curves 
indicates that the cutoff function C increases first with 
increasing argument, before it shows the exponential de- 
cay. This bump is believed to contain all the trees that 
would sit in larger clusters if the system was exactly at 
the critical point | |T5[ |. 

Figure ^ shows the fire size distribution for parameter 
values such that the mean density is ten percent above 
its critical value. As discussed in the previous subsection, 
system spanning fires occur, and their size scales as L^. 
These fires are responsible for the peaks in the fire size 
distribution. Similar peaks occur in equilibrium critical 
systems in the ordered phase, for example in the cluster 
size distribution of a percolation system above the perco- 
lation threshold. In a percolation system, the occurrence 
of such peaks implies that the system size is larger than 
the correlation length, which is identical to the cutoff in 
the radius of the finite (i.e., non system spanning) clus- 
ters. In our system, however, we do not see such an 
exponential cutoff to the size distribution of the finite 
clusters. Instead, the curve for L = 800 in figure ^ ap- 
pears to obey a power law from s ~ 100 up to the point 
where the peak begins. The explanation for this unusual 
behaviour must lie in the large temporal fluctuations in 
the density. The density is only sometimes so large that 
the large fires, which have a size of the order occur. 



4 



At other times, the density values are different and allow 
for a broad range of other fire sizes. 

The transition between critical scaling and L^-scaling 
can be observed when L is varied for fixed sq, as il- 
lustrated in Figure One can see that the shape 
of the curves changes continuously as L is decreased. 
Clearly, because of this change in shape, finite-size ef- 
fects do not manifest themselves in a scaling behaviour 
sn[s) ~ s^'^'^C{sl L^^). It is impossible to generate a scal- 
ing collapse of different curves, even if their density is 
close to the critical density. Furthermore, as mentioned 
above, the cutoff introduced by the finite system size al- 
ways scales as , due to the occurrence of system span- 
ning fires, and not as L^, as expected for conventional 
critical systems. 

For small system sizes, spanning clusters may already 
occur for densities below pc, an effect which is clearly 
visible in the curve for L = 63. The formation of peaks 
due to finite size effects, was also found in JT^ . 

As mentioned earlier, for conventional critical systems 
a system of small size and a small section of a large sys- 
tem are equivalent. We have argued that this is not true 
for the forest-fire model since the mean tree density and 
the time interval between fires are different in the two 
cases. The next two figures show that also the fire size 
distributions are different. In Figure ^the fire size distri- 
bution of a section of a large system and a corresponding 
small system is shown. The scaling parts and the form of 
the bumps near the cutoff are very different. The fire size 
distribution of a small section of a large system is broader 
than that of a small system. The reason is that a sec- 
tion of a large system can contain a boundary between 
a patch of large tree density and a patch of small tree 
density. This boundary can pass through the section in 
different ways, and the number of trees in the dense part 
can take different values. Since large fires only burn the 
dense part, the size distribution of fires becomes broad. 

Figure ^ shows how the fire size distribution changes 
when smaller and smaller sections of a large system are 
evaluated. Comparing to Figure ^, one sees again that 
the fire size distribution of small sections of large systems 
is different from that of small systems. 

To summarize this subsection, the fire size distribu- 
tion in the presence of finite-size effects does not show 
the features of finite-size effects in conventional critical 
systems. We find a continuous change in the shape of the 
fire size distribution and cutoffs that scale as L^, rather 
than conventional finite-size scaling. Furthermore, the 
fire size distribution in small sections of large systems is 
different from small systems. Our results for the fire-size 
distribution confirm the qualitative transition from a pa- 
rameter region unaffected by finite-size effects to a region 
dominated by system spanning fires that we found in the 
previous subsection. 



C. Probability distribution of the density 

Finally, we studied the temporal fluctuations in the 
values of tree density p{t) by measuring how often a given 
value of p occurs within a sufficiently long time series. We 
denote by w{p)dp the probability that the tree density lies 
in the interval between p and p + dp. The quantity w{p) 
is therefore the probability density for the tree density 
p. We measured p always after the trees were refilled, 
i.e., after step (ii). The results show again a qualitative 
change as the correlation length becomes smaller than 
the system size, reflecting the transition from critical to 
percolation-like behaviour. 

In Figures || to ^ w{p) is shown for different values 
for L and sq. For large enough and fixed sq/L'^, the 
mean tree density increases with increasing system size, 
until it reaches its asymptotic value above pc- For fixed 
L, the mean tree density increases with increasing sg. 
Apart from this increase in mean tree density, the fol- 
lowing other trends are observed: (a) As the mean den- 
sity approaches Pc, the curves for w(p) become broader 
(Figure ||) . This is because the patches of homogeneous 
density visible in Figure |l] become larger with increas- 
ing p, leading to larger global density fluctuations, (b) 
As the mean tree density increases above pc, the shape of 
the distribution becomes asymmetric, with the maximum 
moving from p^ to the percolation threshold Pperc 

~ 0.59 

(Figure The reason is that for p > pc the patchy 
structure is replaced by a more homogeneous (percolation 
like) structure, where the largest fires are system span- 
ning and occur for densities above the percolation thresh- 
old. Once the density lies above the percolation threshold 
the probability that a system spanning fire occurs is very 
high. This is why densities much higher than the per- 
colation occur seldom, explaining the rapid decrease of 
w{p) above the percolation threshold, (c) As the system 
size increases for fixed p = sq/L"^, there occur peaks in 
the density distribution which become sharper and more 
numerous for larger L (Figures |l^ and |ll| ) . This can be 
explained by realizing that the difference between finite 
and system spanning fires becomes more pronounced as 
L increases. In the limit L — > oo, finite fires do not affect 
the density at all, while system spanning fires reduce it to 
a small value. Subsequent filling events then increase the 
density to 1 — exp(— p), 1 — exp(— 2p), 1 — exp(— 3p), etc., 
until the density is above the percolation threshold and 
another system spanning fire can occur. These system 
spanning fires do not always occur at the first instance 
where the density is above Pperc, since lightning might 
strike and empty a site. Also, the density immediately 
after a system spanning fire depends slightly on the den- 
sity before the fire. Therefore, the series of density values 
given above, becomes slightly shifted, depending on the 
density just before the last system spanning fire. These 
shifted series of peaks, in turn, give rise to further pos- 
sible density values above Pperc, leading to an additional 
series of peaks, etc. This is the mechanism leading to 
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the fractal peak structure that emerges as L is increased. 
For smaller L, the effect of small fires leads to a larger 
width of the peaks, which can therefore not be resolved 
when they are close together. 

As mentioned in the introduction, the peaks in w(p) 
are due to the fact that lightning can strike the sys- 
tem only between two filling steps. Had we instead per- 
formed our simulations using a small tree growth prob- 
ability p and a small lightning probability /, lightning 
could strike the system between the growth of any two 
trees. However, such a simulation would be very slow. 
In oder to make sure that our choice of the algorithm 
has no other effect on the results apart from the peaks in 
w{p), we performed a test simulation where so is not 
the same for each filling step. For each filling step, 
we chose sq randomly from an exponential distribution 
P{so) — (i^p)^^ exp(— so/(i^p)). Such an exponential 
distribution results when lightning can strike the system 
between the growth of any two trees with the same prob- 
ability. The mean number of trees growing between two 
lightning strokes is now smaller than before. The reason 
is that the majority of filling steps increase the tree num- 
ber by a value smaller than L^p, and that during large 
filling events the tree density becomes high and most of 
the sites chosen for filling are already occupied. The 
mean tree densities evaluated according to Eq. (jj) are 
consequently slightly smaller than before. 

The probability density 'w{p) for the tree density result- 
ing from this modified algorithm is shown in the insets in 
Figures || and As expected, the peaks have vanished, 
while the change in shape from a curve with peak around 
0.4 to a curve with peak near 0.6 due to finite-size effects 
is the same as before. 



IV. CONCLUSION 

In this paper, we have studied finite-size effects in the 
SOC forest-fire model. As these effects become stronger, 
the system rearranges from a structure with patches of 
different densities to a more homogeneous structure with 
large density fluctuations in time. This rearrangement is 
reflected in the structure of the fire size distribution, in 
the mean tree density, and in the temporal density fluc- 
tuations. Qualitatively similar (although quantitatively 
different) rearrangements are observed when smaller and 
smaller sections of a large SOC system are studied. Due 
to these qualitative changes, conventional finite-size scal- 
ing does not hold. Our work thus demonstrates that 
concepts from equilibrium critical phenomena cannot be 
taken over to the study of SOC systems such as the forest- 
fire model. Instead, these nonequilibrium critical systems 
show generically new features unknown in equilibrium. 
As the scaling ansatz eq. ^ which is based on a single 
length scale ^ Sq '^^'^ only approximately correct, 
the true asymptotic scaling behaviour of the model is still 
an open question. 



We suggest that the reason for the unconventional be- 
haviour of the SOC forest-fire model is the fact that two 
qualitatively different types of fires occur: those that 
burn down a patch of high tree density of fractal di- 
mension 2, and those that burn down a tree cluster of 
a smaller fractal dimension within a region of a tree den- 
sity below the percolation threshold. As a consequence, 
the scaling behaviour of the system cannot be charac- 
terized using only one length scale. While the super- 
position of the two types of fires creates the impression 
of simple scaling as long as finite-size effects are small, 
the difference between them becomes clearly visible for 
smaller system sizes, where system spanning fires receive 
a larger weight. We suggest that the superposition of the 
two types of fires is also responsible for the other uncon- 
ventional features of the SOC forest-fire model listed in 
the introduction. 



Models related to the present one have been studied in 
[p|,p5[ and 0. In these models, the tree density is globally 
conserved by filling exactly the same number of trees into 
the system that have been burnt. As long as the density is 
below the critical value, these models are equivalent and 
show the critical behaviour of the SOC forest-fire model 
as the critical density is approached from below. They 
were introduced for the purpose of studying the SOC 
forest-fire model beyond the critical point, i.e. for den- 
sities larger than the critical density pc- As the density 
increases beyond the critical density, both models un- 
dergo large-scale rearrangements. In [^,^, where trees 
are refilled only after the end of a fire, the new structure 
consists of a finite number of large domains of different 
density. In ||^, where each tree is refilled into the sys- 
tem immediately after it is burnt, the new structure has 
a continuously burning fire, and resembles the forest-fire 
model without lightning introduced earlier by Bak, Chen, 
and Tang , which shows spiral-shaped fire fronts . 
In both these models, the dynamics in the restructured 
state are dominated by large fires burning forests of a 
fractal dimension two, similarly to the restructuring due 
to finite-size effects reported in this paper. 



Let us conclude by noting that it is unclear whether the 
behaviour in higher dimensions resembles that in two di- 
mensions. Clearly, as long as the "patchy" structure with 
two qualitatively different types of fire occurs, mean-field 
theory which neglects all spatial structure ||l7|-[l9[| cannot 
apply, and the system must be below its upper critical 
dimension. The recent paper by Broker and Grassberger 
p8[ on the forest-fire model without lightning indicates 
that unusual scaling behaviour can occur also in 3 and 
4 dimensional forest-fire models. If 6 is the upper crit- 
ical dimension of the forest-fire model, as suggested in 
P, p7|p8| , then the scaling behaviour of the forest-fire 
model should be conventional above 6 dimensions. 
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FIG. 1. Snapshot of the SOC forest-fire model for p c± 
40.8% and L=4096. Trees are black and empty sites are white. 
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FIG. 3. Scaling collapse of the fire size distribution for sys- 
tems with the parameters sq/L^ = 0.001 and L = 1600, 800, 
400, 200. The measured mean tree densities are p = 0.40, 
0.395, 0.38, 0.36. C is a suitable scaling constant for each 
curve. 
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FIG. 2. Lines of constant mean tree density in a 
So vs. plane. The bold solid line repre- 

sents the separatrix between the SOC and the percola- 
tion-like behaviour. The remaining lines represent constant 
p = 0.47,0.455,0.43,0.42,0.0403,0.40,0.35,0.343 (from top 
to bottom). The line for p = 0.40 is derived from interpo- 
lated results. 
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FIG. 4. The size distribution of fires for so/i^=0.15 and 
p = 0.454, for the system sizes L = 10, 100, 200, 400, 800 (as 

the peaks move from left to right). The circles mark the points 
of exact scaling, taking the peak of the L = 800 curve as 
reference. 
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FIG. 5. Transition from critical to porcolation-like be- 
haviour as the system size L is decreased for the fixed pa- 
rameter so=200. The parameter L of tlic curves are from 
right to left: L = 1300, 100, 63, 50, 40, 20. The measured 
mean tree densities are p = 0.385, 0.392, 0.402, 0.414, 0.432, 
0.577. 
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FIG. 7. Fire size distribution for systems with L = 800 
(dot-dashed line) and of subsections of this system of size 
I =100 (solid line) and I =20. We used sq/LP' = 0.001 and 
measured the mean tree density p « 0.40. 
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FIG. 6. Comparison of the fire size distribution of a section 

of size / =63 of a larger (L =600, sq =1440) system (solid lino) 
and a small system with L =63 (so =200) (dotted line). For 
both systems we measured a mean tree density of p =0.401. 
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FIG. 8. w(p) for systems with the parameters S()/L^=0.005 
and L = 63, 100, 1600 (as the peaks move from left to right). 
The measured mean tree densities are p = 0.345, 0.367, 0.407. 
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FIG. 9. w{p) for systems with the parameters so/ L^— 0.05 
and L = 63, 1600 (as the peaks move from left to right). The 
measured mean tree densities are p = 0.402 and 0.423. The 
inset shows w(p) for the test simulation with an exponential 
probability distribution for so- The values of L, and the mean 
value of So arc the same as in the main figure, the measured 
mean tree densities are p = 0.39 and 0.408. 
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FIG. 10. for systems with the parameters so/i^ =0.15 
and L = 115 (for the solid line), L = 1600 (for the dotted 

line), L = 10 (for the dash-dotted line) The w{p) values for 
the L = 10 system are multiplied by 2. The measured mean 
tree densities are p — 0.454, 0.455, 0.394. The inset shows 
w(p) for the test simulation with an exponential probability 
distribution for so. The values of L, and the mean value of 
So arc the same as for the solid and dotted curve in the main 
figure, the measured mean tree densities are p = 0.425 and 
0.428. 
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FIG. 11. w{p) for systems with the parameters So/L^=0.20 
and L = 100 (for the sohd line), 1600 (for the dotted line). 
The measured mean tree densities are p = 0.474, 0.473. 
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